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Adiabatically twisting a magnetic moiecuie to generate pure spin currents in graphene 
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The spin orbit effect in graphene is too muted to have any observable significance with respect to its application 
in spintronics. However, graphene technology is too valuable to be rendered impotent to spin transport. In this 
communication we look at the effect of adiabatically twisting a single molecule magnet embedded in a graphene 
monolayer. Surprisingly, we see that pure spin currents (zero charge current) can be generated from the system 
via quantum pumping. In addition we also see spin selective current can also be pumped from the system. The 
pure spin current seen is quite resilient to temperature while disorder has a limited effect. Further the direction of 
these spin pumped currents can be easily and exclusively controlled by the magnetization of the single molecule 
magnet with disorder having no effect on the magnetization control of the pumped spin currents. 


I. INTRODUCTION 

Graphene has been a revolutionary material of this the 21st 
century. Its remarkable that since its discovery around a 
decade ago it has captured the interest of the scientific com¬ 
munity in ways that even High Tc superconductors in its hey 
days couldn’t. A cursory look at graphenes properties would 
provide reasons for this. Graphene is light with huge tensile 
strength and crucially in contrast to most materials, electronic 
energy is linearly proportional to its wave vector and not its 
square and it shows quite remarkable quantum phenomena 
like Klein tunnelingQ) and room temperature quantum Hall 
effect^, etc. Notwithstanding the many advantages of using 
graphene to transport electrons faster and with less dissipa¬ 
tion it has been observed that spin orbit effect in graphene 
is at best negligible thus making them of no practical use for 
applications in spintronics. Then how do we exploit the mani¬ 
fest advantages of graphene and marry them to spintronics. In 
this communication we address this issue. We intend to em¬ 
bed a monolayer of Graphene with a single molecule magnet. 
We then adiabatically modulate two independent parameters 
of the resulting system to pump a pure spin current. The two 
parameters we modulate are one- the Magnetization strength 
of the magnetic molecule which has been shown to be con¬ 
trolled by twisting iti. Secondly, we add a delta function like 
point interaction, which can either be an isolated adatom in the 
Graphene monolayer or an extended line defect or even a very 
thin potential barrier with the strength of the potential barrier 
being controlled by a gate voltage. Using these two parame¬ 
ters we intend to pump a pure spin current. 

Adiabatic quantum pumping is a means of transferring charge 
and/or spin carriers without applying any external perturbation 
(like voltage bias, etc.) by the cyclic variation of two inde¬ 
pendent system parameters. The theory of adiabatic quantum 
pumping was developed in 1 9980 based on the fact that adia¬ 
batic modulations to a 2DEG can lead to changes in lecal den¬ 
sity of states and thus current transport. In 1999, an adiabatic 
quantum electron pump was reported in an open quantum dot 
where the pumping signal was produced in response to the 
cyclic deformation of the confining potential(i|. The variation 
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Figure 1. A single molecule magnet at x = 0 is adiabatically twisted 
along with an electrostatic delta potential at x = a to generate pure 
spin currents in graphene. 


of the det’s shape changes electronic quasi energy states of 
the det thus ‘pumping’ electrens from one reservoir to another. 
The AC voltages applied to the quantum dot in order to change 
the shape result in a DC current when the reservoirs are in 
equilibrium. This non-zero current is only produced if there are 
at least two time varying parameters in the system as a single 
parameter quantum pump dees not transfer any charge. Later 
on the study of adiabatic pumping phenomenon has been ex¬ 
tended to adiabatic spin pumping both in experiment as well as 
theory|l,[3|. In the experiment en quantum spin pumping one 
generates a pure spin current via a quantum dot by applying 
an in plane magnetic field which is adiabatically modulated te 
facilitate a net transfer of spin. Many spintronic devices, such 
as the spin valve and magnetic tunneling junctions^, are as- 
seciated with the flow of spin polarized charge currents. Spin 
polarized currents coexist with charge currents and are gener¬ 
ated when an imbalance between spin up and spin dewn car¬ 
riers is created, for example, by using magnetic materials or 
applying a strong magnetic field or by exploiting spin-orbit cou- 
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pling in semiconductors^. The ^in transport aided by bias 
voitage have been studied eariierig-nlli. 

More recentiy, there has been an increasing interest in the 
generation of pure spin current by quantum pumping without 
an accompanying charge currentlg. In monoiayer graphene, 
a group has attempted to address this issue by propos¬ 
ing a device, consisting of ferromagnetic strips and gate 
eiectrodes[l3- However, ferromagnetic strips in graphene 
are unwieidy and iead to further complication as spin poiar- 
ization is iost at the feromagnet-graphene interface. In this 
work we look at only a piece of monoiayer graphene without 
any ferromagnets but with two impurities one magnetic(SMM) 
the other non-magnetic. Independentiy moduiating these two 
we show generation of pure spin currents. The generation 
of a pure spin-current is oniy possibie if aii spin-up eiectrons 
fiow in one direction and equal amount of spin-down eiectrons 
fiow in the opposite direction, in this case the net charge 
current Icharge = h+h vanishes whiie a finite spin current 
hpin = h — h exists, because = —/;, where or /; are the 
eiectron current with spin-up or spin-down. Further we aiso 
see generation of spin seiective currents nameiy either oniy 
or /; and this impiies the totai magnitude of charge current is 
identical to that of the spin current. 

in this work we show how pure spin currents and spin seiec¬ 
tive currents can be generated via twisting a singie moiecuie 
magnet. Singie moiecuie magnets are of great importance for 
moiecuiar spintronics. Such moiecuies can now be synthe¬ 
sized, for exampie the moiecuie Mn^ac, which has a ground 
state with a large spin (S=10), and a very siow reiaxation 
of magnetization at iow temperate reO^]. SMM’s do not just 
have integer spin values a number of SMM’s like Mn403 
complexllll have large half integer spin S = 9/2. However, in a 
literal twist to this tale (pun intended), a few years back scien¬ 
tists revealed that these big spin molecules can be twisted and 
thus their magnetization can be easily altered El. This fact we 
utilize in our adiabatic pumping mechanism. The mechanism 
of twisting is not difficult to implement. SMM’s have a metallic 
core (for example in Mnn complexes the metallic core con¬ 
sists of theMn-(9—A^-Mn moiety) and a deliberate targeted 
structural distortion of it leads to change in the magnetization 
of the SMM so much so that a initially ferromagnetic SMM can 
undergo a twist to a anti-ferromagnetic SMM. The twist can be 
effectively controlled by two methods-(i) A simple substitution 
of one type of atom (in Mnu complexes a H atom) by another 
sterically demanding atom (again in Mnu complexes by a Me, 
Et, Ph, etc. atom) or by (ii) the other way this can be achieved 
is through the use of hydrostatic pressureEl- In this method 
external hydrostatic pressure is applied on the SMM leading to 
a twisting of the metallic core and thus a net change in mag¬ 
netization. While the substitution method leads to a discrete 
change in the magnetization the pressure method leads to a 
continuous change. Therefore, the pressure method is much 
more amenable to the quantum pumping mechanism wherein 
the pressure applied could be continuously changed albeit adi- 
abatically leading to a adiabatic change in magnetization and 
thus the exchange interaction J. 


II. THEORY 

Graphene is a monatomic layer of graphite with a honeycomb 
lattice structure 01 that can be split into two triangular sublat¬ 
tices A and B. The electronic properties of graphene are effec¬ 
tively described by the massless Dirac equation. The presence 
of isolated Fermi points, K and K', in its spectrum, gives rise 
to two distinctive valleys. We consider a sheet of graphene on 
the.x-y plane. In Fig.[T|we sketch our proposed system. For a 
quantitative analysis we describe our system by the massless 
Dirac equation in presence of an embedded single magnetic 
molecule(SMM). The Hamiltonian used to describe a SMM has 
the following terms: 

Hsmm = -DSI-Js.S ( 1 ) 

The first term represents the Energy of the SMM. D is an uni¬ 
axial anisotropy constant and is the z component of the spin 
of the Molecular magnet. The second term is most relevant 
to us since we deal with electron transport. The Dirac elec¬ 
trons in graphene interact with SMM only via the exchange 
term —7s.S. Further the magnitude in realistic SMM of D 
is very small as compared to 7, D = 0.292Kelvin = 0.025 
meV while J = 100 meV almost a thousand times laroerfll. 
Though anisotropy is very important term in preserving the in¬ 
trinsic properties of SMM, it’s effect on conduction electron is 
hardly a small correction to energy, nothing else. The first term 
will only be relevant in electronic structure calculations for elec¬ 
tronic transport calculations it is not relevant and the only term 
of interest in the exchange coupling. Thus we only considerthe 
second term in the subsequent analysis. We consider a single 
molecule magnet at .x = 0 and another electrostatic delta po¬ 
tential atx = a nm. The Hamiltonian of a Dirac electron moving 
along x-direction can be written as 

H = tiVFO.p + Js.Sd{x)+Vd{x — a). (2) 

The first term represent the kinetic energy term for graphene 
with o the Paulli matrices that operate on the sub-lattices A 
or B and p = ipx,Py) the 2D momentum vector, second term 
is electron interacting with single molecule magnet and third 
term is an electrostatic delta potential. In the second term 
J represents the exchange interaction which depends on the 
magnetization of the SMM and twisting the SMM changes the 
magnetization thus effectively changing J. s represents the 
spin of the Dirac electron while S represents spin of the SMM. 
V represents strength of the adatom situated at a distance a 
from SMM. Here ti^vp (set equal to unity hence forth) are the 
Planck’s constant and the energy independent Fermi velocity 
for graphene. Eq.[2]is valid near the valley K in the Brillouin 
zone and ‘P = [\|/^(r),\|/f (r)] is a spinor containing the elec¬ 
tron fields in each sublattice. The Hamiltonian for /T'-valley can 
be obtained by replacing py by (—py) in Eq. (2j. Because of 
this symmetry i.e; HK{py) = Hi^i{—py), transport coefficients 
will be same in both valleys. So we are confining our discus¬ 
sion only in /T-valley. 

To calculate the quantum spin pumped currents we need to 
introduce the basic theory of quantum pumping which is quite 
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well known as well as solve the scattering problem for electron by If the amplitude of oscillation is small, i.e., for 

with spin (up or down) incident from either left or right. We sufficiently weak pumping {5X,- ^X,), we have, 

introduce them below: 


A. Quantum pumped currents 




ewdXi5X2sin{Q) y 

dX2 


( 8 ) 


To calculate quantum pumped currents we proceed as follows: 
Thus charge passing through lead p- to the left of Molecular 
magnet, due to Infinitesimal change of system parameters is 
given by- 


dQ.,it)=e[^mit) + ^dX2it)] (3) 
with the spin current transported in one period being- 

^ _ ew dXi ^ dNa^dX2^ 

2nJo ^ dXi dt ^ dX2 dt ' ^ ’ 

In the above t = 2%/w is the cyclic period. The quantity 
dNa^/dXi is the emissivity which is determined from the el¬ 
ements of the scattering matrix. In the zero temperature limit 
by- 


dNa^ _ 1 y„, ^^pv go/H,'. 

ax, ^ > 


dXi 


(5) 


Here denote the elements of the scattering matrix as de¬ 
noted above, as evident ,u,v and i can only take values 1,2, 
while a.af takes values t and i depending on whether spin is 
up or down. The symbol “3” represents the Imaginary part of 
the complex quantity inside parenthesis. 

The spin pump we consider is operated by adiabatically twist¬ 
ing the magnetic molecule thus modulating the magnetic inter¬ 
action between the magnetic molecule and scattered electrons 
“J” and the strength of the “delta” function potential V, herein 
Xi = J = Jo + Jpsin{wt) andX 2 =V = Vo + VpSin(wf+ 9). A 
paragraph on the experimental feasibility of the proposed de¬ 
vice Is given above the conclusion. As the pumped current is 
directly proportional to w (the pumping frequency), we can set 
it to be equal to 1 without any loss of generality. 

By using Stoke’s theorem on a two dimensional plane, one can 
change the line Integral of Eq. 0 Into an area Integral, see for 
details Ref|^ 


Ia^, = e [ dXidX2[z— 
Ja oXi 


d dN, 


ail 


d dN, 


ail. 


dX2 dX2 dXi 


( 6 ) 


Substitution of Eq. 0 Into Eq. |6) leads to. 


^aii — C 


dXidX2 


o/=t,^;v=1.2 


^ dXi dX2 


(7) 


This current is for a particular angle of Incidence ((|)), as the 
scattering amplitudes depend on (|). So hereafter, we replace 


In the considered case of a magnetic molecule and delta func¬ 
tion potential the case of very weak pumping is defined by: 
Xp <^Jo = Vo), and Eq. 8 becomes- 


^aiii*^) 


k L 3( 

o/=t.t.;v=1.2 


ds-r 

dJ dv ^ 


(9) 


wherein. 


lo 


ewx^p sin( 0 ) 
2 % 


As we consider only the pumped currents into lead 1(left of 
magnetic molecule), therefore ,u = 1. Further we drop the p 
Index In expressions below. 

For weak pumping, we have the total pumped up-spin current 
given as follows: 




,n. 


,tt-* , 




I q 7 ^12 

V ay dJ J 


ds 


dv dJ 


( 10 ) 


Similarly, we can calculate the spin down current for the case 
of weak pumping by replacing t—5-i and vice-versa. 

For strong pumping, we have the total up-spIn current given 
as: 


ew 




dNA; dXi 
dX\ dt 


dN^dX2, 


with 


dNap 

dXj 


1 
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12 ) 


+ 3 (^ 4 -)]. 


( 11 ) 


here we have dropped the lead index since we pump always 
to left lead (left to SMM) and Is complex conjugate of 
'. To obtain the total current, we integrate over (|). Then the 
total spin-up pumped current for both weak as well as strong 
pumping becomes: 



( 12 ) 


Similarly, we can calculate the spin down current for the case 
of weak pumping by replacing f—5^4- arid vice-versa. In the 
above equations the scattering amplitudes represent- 
reflection amplitude when spin-up electron Is coming from the 
left side and reflected to the spIn-up state. 
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ijj = reflection amplitude when spin-down electron is com¬ 
ing from the left side and reflected to the spin-up state 

= f-j., transmission amplitude when spin-up electron is com¬ 
ing from the right side and transmitted to the spin-up state. 

= f|, transmission amplitude when spin-down electron is 
coming from the right side and transmitted to the spin-up state. 
Numerically, we have calculated and f| and substi¬ 

tuted in the above expression to obtain the spin-up pumping 
current. 

Effect of finite temperature: So far our discussion is for zero 
temperature. The temperature effects could be easily captured 
by multiplying a factor [-df{E)/dE] with pumping current and 
integrating over energy of incident electron[T3 i.e; 

Ig= [ ^ [ [--§^]laWcos{ft))dEdft),a=t,i. {t3) 
J-n/lJO dE 

where f{E) is the Fermi-Dirac distribution function. 


B. Solving the scattering problem 

Let us consider an incident spin- up electron from left of 
magnetic molecule {x < 0) with energy E. The electron can 
be reflected or transmitted to spin-up/down state. Then the 
spinors, for the angle of incidence (|), in the various regions are 
given as: 

Region-1 {x < 0): 


¥aW 



r^e '*^)x 


m 




-rie-‘^^-‘hn,+i 


Region-ll (0 < .x < a): 






(14) 


(15) 


(16) 


¥bW 




and in Region-Ill (x > a): 


(17) 


<(x) = 




(18) 


¥f(x) 




(19) 


k = EfCos{(^) with Ef{> 0) being the Fermi energy, x’s 
denote the eigen states of the z-component of spin 
operator for SMM, S^Xm = rnXm with m being the corre¬ 
sponding eigen value. The scattering is elastic and the 
exchange interaction conserves the z-component of the 
total spin 5-Fi. The exchange operator in the Flamiltonian, 
s.S = (l/2)(i^5+ -Fi+5^) acts as a spin-flipper for 

electrons with different values of to those of while for 


same values it acts as a normal barrier. 


S-S+ 


1 

0 


Xm — 


-[?] 

Xm+\ and s+S 

0 

1 

X 

3 

II 

1 

0 

Xm-1 With 

F = \/{S — m) (5 -F OT -F 1) 

and F' = yJ{S + m 

){S — m+ 1). 


Flere, ± Sy and = Sx± Sy are the raising and 

lowering operators. 


In solving the scattering problem from a delta potential the fol¬ 
lowing two boundary conditions have to be met: i) continuity 
of the wave functions at the boundary and ii) discontinuity of 
the energy at the boundary. Flowever, these boundary condi¬ 
tions work only if the system is described by the Schroedinger 
equation, but does not work for Dirac material like graphene. 
For Dirac equation, the wave functions on either side of delta 
potential are not continuous across the boundary. 

At x=0, while taking integration on the both sides of the Dirac 
equations, H\\f = Eij/, one is stuck at the following integration: 

¥(0) = / i|/8(x)c/x (20) 

Jx=o- 

because of the discontinuity of wave functions at the boundary. 
Flowever, there is a standard way to avoid this difficulty, which 
has been widely used by many authors, by taking the average, 
i.e.. 


¥(0) = ^[¥(jc = 0+)-f\(/(x = 0 )], (21) 

where the delta function potential mentioned is not an exact 
delta function potential but a point like interaction lioHiill . The 
above idea can be deployed at the boundary x=0, which leads 
to the two equations as 

- ihvFiVI= 0^) -Vb{x = 0^)] 

= '^^■S[VI {x = 0+) +i(/^(x = 0^)j22) 

- i^vfIVI= 0^) -Va{x = 0^)] 

= {x = 0+) +\(/b(x = 0^)j23) 

After substituting the wave functions in Eq. j22) and Eq. j2^, 

a^{iam + e‘'^) +a^{iaF) +b^{iam — +b^{iaF) 

~F r-|-(/aOT -F + r^{iaF) = e'^ — iam, (24) 

a^{iaF) + a^[e“‘^ — ia{m+ 1)] +b'f{iaF) —b^[ia{m+ 1) -Fe^'*^] 
-F r^{iaF) + — ia{m+ 1)] = —i(XF, (25) 


In the above equations, t arid 4- stand for spin-up and 
spin-down electron. Flere, and f|,f; are the re¬ 

flection and transmission amplitudes respectively. Also, 
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(a) Quantum pumping current Vs. Energy of incident electron for weak 
pumping 


(b) Quantum pumping current Vs. strength of the delta like point interaction for 
weak pumping 




(o) Weak Pumping: Quantum pumping current Vs. Strength of the molecular (d) Strong Pumping: Quantum pumping current Vs. Strength of the molecular 

magnet magnet 


Figure 2. Pumped spin currents at zero temperature, other parameters are mentioned in Figure. 


a^{l+iame‘'^)+ a^^{iaFe‘'^)+b^{l—lame —b^{iaFe '‘I’) 
— r|(l + iame^’'^) — r^{iaFe^''^) = 1 — iame'^, (26) 

and 

a^{ia.FF^) +a^[\ —ia{m+ 1)6'*^] 

— b^{iaFe^"^) +ia{m+ 

-r^{iaFe-'^)-ri[l-ia{m+\)e-'^] = iaFe'^. (27) 

Here, a = J/{2 Hvf). The mixing of the spin-up and spin-down 
components in the above equations are attributed to the ex¬ 
change operator s-S. 

At x=a, the boundary conditions are given as: 

— r/iVf [\|/™(x = fl+) — ij/^ (x = a^)] 

= y = a+) -f i|/^^(x = a^)] 


and 

— r/tvf [\|/™(x = fl^) — i|/^^ (x = a“)] 

= '^WB{x = a+)+^i^^i{x = a-% 

iead to 

{iG - e'^) + {iG + e-‘^) + + iG) = 0, 

(28) 

{iG - e'^) + {iG + e-'^) + {F^ + iG) = 0, 

(29) 
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-F/Ge'*'’) - +/Ge-'‘f) +fte''^“(l +iGe‘‘f) = 0, 

(30) 
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and 


aie‘^%-l+iGe‘^)-b^e-‘’‘‘‘{l+iGe-‘^)+t^e‘’‘‘‘{l+iGe‘^) = 0 , 

(31) 


where, G = V j{Itivp)- 

The eqns. 1241311 contain 8 unknown probability amplitudes, 
which can be solved numerically to confirm | ft 1^ + I ''t P + I 
P + I P= 1- Similarly for the case of spin-down inci¬ 
dent electron from the left side, we can get scattering ampli¬ 
tudes. This procedure can be repeated appropriately for spin- 
up/down electron coming from right side. 


h.S 


aT 


ikx 


Bje-*''" 

B i ^-ikx 
0 = 



t Jkx 



I o-iVcx 


pi p-ikx 


Figure 3. A typical diagram of the disorder added device in which 
SMM and adatom is separated by a distance a nm. 


III. RESULTS AND DISCUSSIONS 


Different components of quantum pumped currents i.e; spin- 
up (/t), spin-down (/;) and net spin current (h = — and 

charge current {G =/t + 4) obtained numerically using 
Eq. 12 at zero temperature as shown in figures 2 (a)-(d). 
Figure 2(a) and 2(b) show the quantum pumped currents 
versus energy of the incident electron (£) and 
strength of adatom/line defect for weak pumping, without dis¬ 
order. The following numerical parameters are used (men¬ 
tioned in the figure also): the spatial separation between SMM 
and the adatom a = 10 nm, the strength of exchange inter¬ 
action between electron and the molecular magnet J = 0.3 
eV-nm, m = — 1 /2 and spin of the molecular magnet S = 3/2. 
The strength of line defect potential V = 0.3 eV-nm and the 
strength of the time varying modulations on top of V and J are 
taken as Xp = 0.05 eV-nm in strong pumping case. From fig¬ 
ure 2(a), it is seen that for a suitable energy E (= 54 meV), the 
total charge current completely disappears leaving behind a 
pure spin-current. Fig. 2(b) shows that individual spin pumped 
currents slowly varying with the increase of ordinary potential. 
The is decreasing while /; is increasing with V. The pumped 
currents in weak pumping regime are in units of /q as in Eq. 8, 
while in strong pumping regime are in units of ew/27t. Another 
important point in the weak pumping case is that, by suitably 
choosing E, I/, can be completely suppressed to pump only 1 1 
a spin up selective current. 

Figures 2(c) and 2(d) are plotted to check how the strength of 
the exchange interaction (J) affects the different pumped cur¬ 
rents. We have taken energy of incident electron {E = 54 meV) 
for which pure spin current was observed. From figure 2(c), 
we see that the pumped currents attain a maximum around 
J = 0.1 and —0.13 eV-nm and then start diminishing towards 
zero, and pure spin current appears at 7 = 0.3 eV-nm, and 
spin selective current /; appears for J > 0.5 eV-nm. However, 
in strong pumping as shown in figure (d), pure spin current ap¬ 
pears at the same position at 7 = 0.3 eV-nm. For both weak 
and strong pumping plots, we have chosen the phase differ¬ 
ence between two modulations as %/2. One important fact 
which is clearly noticeable is the fact that ‘J’ acts as a current 
switch. By changing J from positive to negative, all the pumped 
currents change sign. This shows that the magnetization of 
SMM controls the direction of current flow in graphene. Since, 
we aim to control this by the twist, effectively twisting the SMM 
changes the direction of spin currents. This is a key result of 


our work. 


A. Effect of Disorder and temperature 

In this subsection, we discuss the effects of disorder distributed 
randomly in the system. We have modeled the present device 
in such way that SMM and adatom/line defect are at the ex¬ 
treme ends of the sample and disorder is confined between 
SMM and adatom/line defect. Here each disorder is consid¬ 
ered to be a delta potential (point interaction as mentioned be¬ 
fore). We solve this problem by using transfer matrix approach. 
The presence of multiple delta potentials causes a number of 
confined regions between SMM and the adatom. The general 
form of the wave function in each region can be written as: 


i|/„(x) 




(32) 


Here n = 0,1,2,3...(A^— 1),A^ corresponding to different re¬ 
gions bounded by the delta potentials, as shown in Fig. 3. 
The above wave function is for A-sublattice, the phase factor 
exp(±/(|)) is multiplied with transmission (reflection) amplitude 
to get the same for B-sublattice. The next step is to find out 
the total transfer matrix which connects the wave function am¬ 
plitudes between extreme left and right. To do so, first we find 
transfer matrix across SMM i.e; between region ‘n = 0' and 
= 1' as in Fig. 3, 


A 

A 

B 

B 


t 

1 

1 

1 

t 

1 

; 

1J 




r a’^ 1 
^0 

dI 

X>r\ 


(33) 


where mI*’**! is the transfer matrix across SMM, expressed as 
=Mi/Mo with 


Mo = 


^-igm -igF igm-^c 
-igF ^ + ig{m+l) -igF 

1 - igm^ -igFi, 1 + igm'i.c 
-igF^ l+ig{m+l igF^c 


-igF 

-ig{m+l)-^c 

igF^c 

l-ig{m+l)^c . 
(34) 
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(a) Quantum pumped currents in the weak pumping regime vs. spin of the (b) Quantum pumped currents in the weak pumping regime vs. J for 

SMM for disordered graphene. disordered graphene. 

Figure 4. Pumped spin currents in presence of disorder. 



0/2JI 


Figure 5. Quantum pumping current Vs. phase angle between two 
modulations for strong pumping. Solid line is for zero temperature 
case while dashed line Is for finite temperature. 


transfer matrix between ‘n = N' and ‘« = At — 1' as 


A 

A 

B 

B 


t ■ 
N 
; 

N 

t 

N 

4 

N . 




A 

A 

B 

B 


t 

N-\ 

N-\ 

t 

W-1 

N-l 


(36) 


where '1 is the transfer matrix across adatom, ex¬ 
pressed = Mn/Mn^\ with 


Mn-i 


^-iG 0 
0 ^ - tG 

1-iG^ 0 

0 1 - /G^ 


iG-^c 0 1 

0 IG-^c 
l+IG^c 0 
0 1+iG^c. 


(37) 


and 


Mn = 


^ + iG 0 IG-^c 0 
0 ^ + iG 0 IG-^c 

l+IG^ 0 1-iG^c 0 
0 1+tG^ 0 l-iG^c 


(38) 


and 



^ + igm 

igF 

igm - 

igF 

]\4^ _ 

igF 

'i,-ig[mF\) 

igF 

-ig{m+\)-'i,c 

IVI \ — 

1 + igm^ 

igF^ 

1 - igm'i.c 

-igF^c 


- igF^ 

l-ig{m+l)^ 

-igF^c 

l+ig{m+l)^c 





(35) 

with ^ 

— exp(i(|)) and = exp(—; 

(|)). Similarly we can get the 


transfer matrix across any ordinary potential, for example, the 


Since the adatom is modeled as a delta function potential and 
disorder is modeled too as randomly distributed sequence of 
delta potentials with random strengths the transfer matrix for 
any arbitrary interface between adatom and SMM has also the 
same matrix elements as . 

After some straight forward algebra, the connection between 
the wave function amplitudes of extreme left and right is found 
to beliol 


r A^ 1 


r a'^' 1 

< 

= M 

^0 

pi 

L J 


pi 

L “o J 
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A Effect of Disorder and temperature 


IV PUMPING VS. RECTIFICATION 


where 


M = ^ .’“I (40) 


with the propagation matrix between any two succes¬ 

sive delta potential, is given by 


M'jree = 


eikdn 0 0 0 

0 e'“'' 0 0 

0 0 0 
0 0 0 e^'“" 


(41) 


Here, dn is the separation between two consecutive delta po¬ 
tentials. To calculate the reflection and transmission ampli¬ 
tudes, we shall use the scattering matrix (S-matrix), which is 
connected to transfer matrix asfiai 


S = 


1 

(M22)2x2 


[M2\)2x2 1 

det(M) (Mi2)2x2 


M = 


Mil Mi2 

M 21 M 22 


mil mi2 m\2 mi4 
m2\ m22 OT23 "224 
^31 ^32 ^33 ^34 
_ m^i m42 m43 ^44 _ 


(42) 


The reflection amplitude (to the left, as we are calculating 
pumped currents in the left lead) 


^21 ^ 

M 22 '"U 


and transmission amplitude from right to left is 


tr = 


1 

M 22 


fn Hi 

.Ht Hi 


(43) 


(44) 


The tr and r, can be directly used in Eq. GD to obtain the 
quantum pumped currents. We must mention here that disor¬ 
der free pumping current can also be recovered from here by 
using transfer matrix as M = MP’^lM/reeMl'’’^!, where mP’' 1 
and would become the transfer matrix across adatom 
and SMM, respectively. 

Disorder effect is shown in Fig. 4(a). Herein we plot the 
pumped currents as function of the different spin states of 
SMM in the weak pumping regime. There is pure spin cur¬ 
rent for 5 = 3/2, in the figure m = —1/2. To include disorder, 
we have chosen the random spacing between any two poten¬ 
tials in the range 1 — 1.5 nm. The strength of the potential is 
also random and ranges from 50 — 100 meV-nm (dashed line) 
to 100 — 200 meV-nm (dashed-dot line). We have considered 
the pumped currents averaged over 1000 realizations. One 
can see that disorder has a limited effect effect on the pure 
spin current. The position of pure spin current is shifted from 
5 = 3/2 (without disorder) to 5 = 5/2 (disorder:50 — lOOmeV- 
nm) and finally 5 = 7/2 (disorder:100 —200 meV-nm) however 
pure spin currents aren’t killed off. In Fig. 4b we see the ef¬ 
fect of disorder on pumped currents plotted as function of the 
magnetization of SMM. We again see position of occurrence of 
pure spin current (4 = 0) changes from J = 0.3 to 0.45 eV-nm 
as one increases disorder. However, disorder has no effect on 


magnetization switching of pumped spin currents showing the 
resilience of the magnetization switching to disorder. 

Finally, in figure 5 we plot the pumped currents as function of 
the phase difference between modulated parameters. We see 
that pumped current attains maximum at 0 = k/2 and minima 
around 0 = 0,7t. The pure spin current is maintained through¬ 
out the whole range of 0. The temperature effect is shown in 
the same figure, which shows a small damping in amplitudes 
of the individual spin currents. One can see that temperature 
has no noticeable effect on the pure spin currents apart from a 
slight diminishing of the magnitude. To conclude this section, 
pure spin currents in graphene are immune to any temperature 
increase apart from decrease in magnitude while disorder has 
a small effect as it shift the parameter regime for occurrence of 
pure spin currents although it cannot kill it off. 


IV. PUMPING VS. RECTIFICATION 


A major issue which was flagged right from the early days of 
quantum pumping was whether the Switkes experiment!^ was 
a real demonstration of quantum pumping, since the pumped 
current was observed to be symmetric with respect to Mag- 
netic field reversal just like the two terminal conductance!^. 
However since pumped currents are functions of scattering 
amplitudes and not scattering probabilities they should have 
no particular symmetry with respect to magnetic field reversal 
unless the system itself had some particular svmmetrvli^. As 
the Switkes expt. system did not possess any particular sym¬ 
metry it was quickly recognized that the current attributed as 
a a pumping current was in effect a rectified current which de¬ 
pend on the Conductance of the svstemli^l. However there 
could have been a pumped current which was masked by the 
rectified currents. The origin of rectified currents is because 
experimentally at the nanoscale it is difficult to control time 
varying parameters. Most naturally time varying parameters 
couple to input and output leads and instead of only pumping 
a current there is in addition a transport current defined the 
net conductance through the system. So any quantum pump¬ 
ing at the nanoscale will have rectified currents and there¬ 
fore it become imperative to have a scheme to differentiate 
these currents. The rectified spin up current is defined as: 

Irect = IsdXidX 2 {Ci^ — C 2 ^),Xi’s are the modulated 
parameters. In the weak pumping regime we have- Irect = 
Io{Ci^-C 2 ^), with /o = we^sm{Q)dXidX2R/4K^n. 

In Fig. [6]We can see clearly that the Conductance (both spin- 
up as well as spin-down) are symmetric with respect to small 
values of J. Thus unlike pumped currents whose direction can 
be changed by changing the magnetization from positive to 
negative, the conductance on the other hand shows no such 
effect. The pumped currents are completely asymmetric as 
function of J as seen in Fig. 2(c) and (d). Thus even if rectified 
currents will be present in the system the pumped spin current 
will be distinguished because of their asymmetric nature with 
respect to magnetization reversal. 
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Figure 6. Conductance (spin up and spin down) vs. J. For smaii J 
values G’s are symmetric. Compare with Fig 2(c) and (d). Pumped 
currents up and down spin are asymmetric. 


V. EXPERIMENTAL REALIZATION AND CONCLUSIONS 

The experimental realization of our pure spin current pumping 
device shouldn’t be too difficult. As already outlined in the last 
paragraph of the introduction of this paper adiabatically mod¬ 
ulating the pressure applied on the single molecule magnet 
would entail a corresponding adiabatically modulated magne¬ 


tization of SMM. The second adiabatically modulated parame¬ 
ter of the device is an adatom placed “a” distance apart from 
SMM. The adatom is modeled as a delta function like point in¬ 
teraction similarly embedded in graphene. A gate voltage ap¬ 
plied to the adatom can change the potential felt by electrons 
scattered from it. When the gate voltage itself is adiabatically 
modulated in time we have all the ingredients for the quantum 
pumping of pure spin currents and spin selective currents. Al¬ 
ternatively, an extended line defect can be created instead of 
adatom, which can be controlled exDerimentallv l23 - [23 . More¬ 
over, one can also use a thin potential barrier experimentally 
which can be theoretically modeled as a delta like potential. 
Similarly, the single molecule magnet is infact a large molecule 
with host of atoms ranging from 30-100 atoms, these atoms 
are arranged not only having vertical but also horizontal ex¬ 
tension i.e., a single molecule magnet will have sufficient ex¬ 
tension along transverse direction. Mention may be made of 
Refliolon single molecule magnets which exemplifies the situ¬ 
ation envisaged. 

To conclude we have shown that notwithstanding the fact that 
spin transport via spin orbit effect is almost impossible to be 
observed in graphene, we have in a novel manner pumped 
pure spin currents and spin selective currents in graphene via 
embedding it with a single molecule magnet. The study of pure 
spin currents in graphene via embedded SMM will be extended 
to spin correlations and whether one can generate entangled 
spin currents which will have potential impact on quantum in¬ 
formation processingfs^ in a subsequent work. 
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